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INTRODUCTION 


Extended time observations (ETO) were made at San Nicolas Island (SNI) 
from March to October, 1987. A small ground station was installed at the NRL 
trailer (Site B) . Hourly averages of air temperature, relative humidity, wind 
speed and direction, solar irradiance, and downward longwave irradiance were 
recorded. The radiation sensors were standard Eppley pyranometers (shortwave) 
and pyrgeometers (longwave). Data records of SNI rawinsonde launches tor this 
period have been requested from the Pacific Missile Test Center (PMTC) . 

This data will be processed in a variety of was to deduce properties of 
the str atocumul us covered marine boundary layer (MBL). From the temperature 
and humidity the lifting condensation level, which is estimate of the 
height of the cloud bottom, can be computed. Combinations of shortwave and 
longwave irradiance statistics can be used to estimate fractional cloud cover 
parameters. 

This paper describes an analysis technique used to estimate the 
integrated liquid water content <W) from the measured solar irradiance. The 
cloud transmittance is computed by dividing the irradiance measured at some 
time by a clear sky value obtained at the same time on a cloudless day. From 
the transmittance and the zenith angle, the cloud LWC is computed using the 
radiative transfer par ameter i z at i ons of Stephens et al. (1984). The results 
are compared with 17 days of mm-wave (20.6 and 31.65 BHz) radiometer 
measurements made during the FIRE IPO in July of 1987. 

CLOUD RADIATIVE TRANSFER COMPUTATIONS 

Stephens (1978) has developed a simple parameterization of cloud albedo 
(Re), transmittance (Tr), and absorption (Ab) in terms of the cloud optical 
thickness (T) and the zenith angle (0). The parameterization is based on the 
two-stream approximation of Coakley and Chylek (1975). A lock-up table of 
values for single scattering albedo (u) and back scatter fraction ( R) is given 
as a function of T and p-cosO). The values are based on results of 
sophisticated 10-layer, 15-band radiative transfer computations. The revised 
look-up table (Stephens et al . , 19B4) was used for the work presented here. 
Non absorbing < X < 0 . 7 pm) and absorbing < X > 0 . 7 pm) wavelength regions are 
considered. The two-stream model is as follows: 

(i) Non absorbing, 4)=1 

Re = S T / X / < I * 0 T/X ) (la) 
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Tr = 1 - Re 


(lb) 


(i i ) Absorbing , w< 1 

Re = <u 2 - 1 ) [exp (Tef ) -exp (-Tef ) 3/R (2a) 

Tr = 4u/R (2b) 

Ab = 1 - Re - Tr (2c) 

where 


u 2 - (1- ui + 28<u)/(l - id) (3a) 

Tef = SQRTC < 1 - u)*(l - id + 2&»)]T/p (3b) 

R = (a + l) 2 exp(Tef) - (u - l) 2 exp(-Tef) (3c) 

The optical thickness is calculated from the cloud W using a parameterization 
obtained from a number of measured cloud properties. 

log lo (T) = 0.2633 + 1.7095 1 og .[ 1 og , „ < W) 1 , X<0.7 (4a) 

1 og i o ( T) = 0. 3492 + 1.6518 1 og.Cl og lo (W)l , X>0.7 (4b) 


Since the Stephens parameterization was developed to apply to some 15 
cloud types, we decided to compare it to a model with microphysics specific to 
the marine str atocumul us , which we shall refer to as the DF6S (Davidson et 
a 1 . , 1984) model. fhe DFBS model is a 15-band, del ta-Eddi ngton radiative 
transfer model that assumes a linear liquid water profile in the cloud with a 
log-radius gaussian cloud droplet distribution where the mode radius is 
determined from the liquid water assuming a constant number of cloud droplets 
(100/cc) . Figure 1 shows a comparison of transmi ttances for the two models at 
three different zenith angles for a generic California stratocumulus with 
cloudtop at 600 m and total water (liquid plus vapor) mixing ratio of 8 g/kg. 
Optical thickness was varied by varying the cloud thickness from 50 m to 450 
m. The comparison is very good but the Stephens par amet er i z at l ons (Eqs. 4) 
appear to underestimate the optical thickness at low values of W, which are at 
the lower limit of his fit. Since marine stratocumulus often have low values 
of W, we used the parameterization given below for W<100 g/m''2 (see Fig. 2), 
which is obtained from the DFBS model: 

logio(T) = -0.35 + 0.91 log lo (W) (5) 


Using this slightly modified form of the Stephens model, a nomogram can be 
constructed in terms of contours of transmittance versus zenith angle and 
integrated liquid water path (Fig. 3). 


RESULTS 

Using the clear sky irradiance, the cloud transmittance was computed for 
each day from July 1 to July 18 for the times 1600 to 0100 GMT. Low incidence 
angles (p< . 4) were not used because of the potential for large errors. For 
the same reason, transmi ttances greater than 0.9 were not considered. For a 
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given value of Tr and p, the value of W was iterated until the Stephens model 
gave the measured transmittance. This computed value of W is compared with 
the values from the NOAA radiometer (Fig. 4). 

The mm-wave radiometer values of W are on average about 507. greater than 
the values estimated from the transmittance. The rms scatter is about 357.. 

The NOAA radiometer is believed to be accurate to 207., but a 507. fractional 
bias is not out of the question. Some of the rms scatter is due to 
differences in sampling geometry, since the pyranometer has a hemi spher i cal 
field of view and the radiometer has a narrow field of view. These 
differences in field of view are somewhat moderated by averaging for one hour. 
Another source of error is the use of a single clear sky irradiance for the 
entire 17 day period, which ignores variability in the meteor ol ogi cal 
conditions above the mixed-layer. Once these issues are sorted out, He expect 
to use this method to analyze the entire data base. 
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Fig. 1 Comparison of transmittance versus integrated cloud INC at zenith 
angles of 0, 60, and 80 degrees for the Stephens (circles) and DFBS (line) 
models. 

Fig. 2 Optical thickness versus integrated cloud LWC for the Stephens 
parameterization of Eq. 4a (solid line), the DFBS calculations (circles), 
Eq. S (dotted line). 

Fig. 3 Contours of transmittance as a function of solar zenith angle and 
integrated cloud LWC from the modified Stephens parameterization. 

Fig. 4 Integrated cloud LUC: the vertical axis is computed from solar 
irradiance versus W from the NQAA radiometer measurements. 
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